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ABSTRACT 

I argue that cool giant companions in most symbiotic binary systems possess magnetic activity 
on a much higher level than isolated, or in wide binary systems, cool giants. Based on the behavior 
of main sequence stars, I assume that magnetic activity and X-ray luminosity increase with rotation 
velocity. I then show that the cool companions in symbiotic systems are likely to rotate much faster 
than isolated, or in wide binary systems, cool giants. The magnetic activity of the cool giant may 
be observed as a global axisymmetrical mass loss geometry from the cool giant (before the hot 
companion influences the outflow), a stochastic mass loss process, i.e., variation of mass loss rate 
with time and location on the giant's surface, and in relatively strong X-ray emission. The variation 
in the mass loss process from the cool giant may cause variation in the properties of jets blown 
by the hot compact companion. I conclude that symbiotic systems should be high-priority X-ray 
targets. 

Subject headings: binaries: symbiotic — circumstellar matter — stars: AGB — stars: magnetic 
fields 
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1. INTRODUCTION 

The level and the nature of the magnetic activity in red giant branch (RGB) and asymptotic 
giant branch (AGB) stars is poorly known and understood. Indications of magnetic activity come 
from polarized maser emission and X-ray emission from a number of AGB stars. Polarized SiO 
maser emission is observed close, typically less than a few stellar radii, to some evolved stars (e.g. 
Kemball &; Diamond 1997), whereas polarized OH maser emission is observed at ~ 10^^ — 10^^ cm 
from the star (e.g., Zijlstra et al. 1989; Szymczak, Cohen, k. Richards 1999; Miranda et al. 2001). 
X-ray emission found in several evolved stars (e.g., Hiinsch et al. 1998; Schroder, Hiinsch & Schmitt 
1998; Hiinsch 2001) may indicate the presence of a hot corona that results from magnetic activity. 

The dynamo in evolved stars, e.g., AGB stars, is different in many respects from the dynamo 
mechanism in the Sun. The main difference is the ratio between the rotation period, Pj-ot = 27r/a;, 
and the convective overturn time Tc (Soker & Zoabi 2002). This is defined as the Rossby number 
Ro = {u)Tc)~^- In the sun and other magnetically active main sequence stars, Ro < 1, whereas in 
evolved stars Ro 3> 1. In cool giant stars, therefore, the amplification is most likely via a turbulent 
dynamo, i.e., an aa dynamo, but where the rotation plays a positive role (Soker 2000b; Soker & 
Zoabi 2002). This mechanism is termed a^u dynamo. If indeed rotation plays a significant role in 
the amplification of the magnetic field in cool giant stars, then systems where these stars rotate 
fast would be expected to be more active. Such systems are symbiotic stars (SSs), which are binary 
systems composed of a mass losing evolved star, basically an RGB or an AGB star, and an accreting 
compact star, in most cases a white dwarf (WD), and in the minority of the systems a close main 
sequence (MS) star (e.g., Belczynski, et al. 2000). 

In §2 I show that the cool companion in symbiotic systems are likely to rotate much faster than 
single cool giant stars, hence they are prime candidates for possessing (relatively) strong magnetic 
activity. The presence of magnetic fields coronae in SSs were proposed in the past, e.g., to account 

for UV emission (Stencel & Sahade 1980). A thorough review of binary evolution and interaction in 
symbiotic stars is given by Iben & Tutukov (1996). I will only concentrate on the expected rotation 
speed of AGB cool companions in these systems, and the expected magnetic activity. In §3 I discuss 
the manifestation of the expected magnetic activity, and in §4 I present my main conclusions, and 
propose some future observations. 



2. SPINNING-UP THE COOL COMPANION 

2.1. Accretion during the main sequence phase 

In this subsection I consider SSs with a WD hot companion. I show that the presently cool 
giant companion was spun-up when it was an MS star and the presently WD was an AGB star. In 
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most SSs the accretion rate into the WD companion is Mace > 10~^Mq yr"^, the orbital period is 
Porb < 10^ day, the mass loss rate (defined positively) from the cool companion is < W~^Mq yr, 
and the WD mass is MwD ^ O.GMq. The parameters of most systems, in particular the mass of 
the WD companion, indicate that the progenitor of the WD companions was relatively massive 
Mp ^ 3Mq. Such stars lose ~ 1 — 2Mq as a slow wind during their AGB phase. Some fraction 
of this mass is accreted by the MS progenitor of the presently giant star. To find the accreted 
angular momentum by the MS progenitor of the present cool giant, I use the Bondi-Hoyle mass 
accretion rate, and the angular momentum entering the Bondi-Hoyle accretion cylinder (Wang 
1981), i.e., having impact parameter b < Ra = 2GMg/v1^ where Ra is the accretion radius, Mg is 
the mass of the accreting star, and Vr the relative velocity of the star and the wind. The accreted 
angular momentum onto the MS progenitor of the present giant is given by (for R^, < a; see relevant 
equations and more references in Soker 2001) 

Jacc = I {^y^MMg + M,)aM^]V2, (1) 

where = MgMp/ {Mg + Mp) is the reduced mass of the binary system, Mp is the mass of the 
WD's progenitor, r/ is the ratio of the accreted angular momentum to that entering the Bondi-Hoyle 
accretion cylinder, and a is the orbital separation (for simplicity a circular orbit is assumed). The 
last term in equation (1) (the []^/^ term) is the total orbital angular momentum of the binary 
system. I took the wind velocity v^, to be equal to the relative velocity Vr = (v^ -|- Vo)^^"^, where 
Vo is the orbital velocity of the accreting star relative to the mass losing star. This is justified for 
the parameters used here and the other uncertainties, e.g., in the value of rj. Scaling the different 
variables with typical values, I find the accreted angular momentum to be 

7 ^107 (v\( M, + MA y^f Mg Y( \f a \ V 

J^,-WJe^^ ^J^ 4Mq J [i.OMqJ [iMRq)[50A\j) Us km s"! j ' ^ ^ 

where J© = 1.7 x lO''^ g cm^ s~^ is the present angular momentum of the Sun. Stars more massive 
than the sun rotate faster, up to 100 times faster, but accrete much more angular momentum, e.g., 
by a factor of ~ 100 for Mg = 3Mq. Of course, the MS can't rotate faster than the break-up 
velocity. The conclusion from the last equation is that the present giant in most SSs reach the 
RGB and AGB with typical angular momentum, hence angular velocity, much larger than isolated 
stars do. 

As an example I consider Mira (o Ceti) a non-SS binary system with an AGB star and a WD 
in a projected orbital separation of ~ 70 AU (Karovska et al. 1997). The true orbital separation is 
larger, but it was smaller before mass was lost by the two components. Taking therefore a = 70 AU 
in equation (2) with all other parameters held the same, gives Jacc — 3Jq. If, on the other hand, the 
slow wind speed were lower, as is quite plausible for a large AGB star, as the progenitor of the WD 
was, say Vu, = 10 km s~-^, then Jacc — 7b Jq. I suggest, therefore, that Mira has reached the AGB 
with an angular momentum larger than typical isolated AGB stars do, despite the orbital separation 
that is too large to allow substantial tidal spin- up (see next subsection). Such an extra spin may 
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enhance magnetic activity (see §3), which may lead to a stochastic mass loss process (§3.2), which 
may explain Mira's asymmetrical shape observed by Karovska et al. (1997) and Marengo et al. 
(2001). 



2.2. Tidal interaction 



The hot companion spins-up the cool giant companion via tidal interaction, whereas mass loss 
from the giant spins it down. I use the equilibrium tidal interaction for giants (Zahn 1977; 1989; 
Verbunt & Phinney 1995) with time scales from Soker (1998). The circularization time of the orbit 
is given by 

r • = 1 2 X 10«1 (^^Y'^' f^^Y' (M^y (M^Y' 
""^^ ■ /, V2OOOL0; \200RqJ \ Mg J \MqJ 

where Lg, Rg and Mg are the luminosity, radius, and total mass of the giant star, M^nv is the giant's 
envelope mass, M2 is the mass of companion that exerts the tide, a is the semi-major axis, and 
/s(e) is a steeply increasing function of eccentricity e (Hut 1982) with /s(0) ~ 1. The transferring 
rate of orbital angular momentum to giant spin is 

<A)rb — ^7~<^orb — j (4) 

ZCL '7~cir 

where Jorb is the total orbital angular momentum of the two stars, and dot stands for time derivative. 
The expression is accurate for a circular orbit; for eccentric orbits there is a factor of correction 
(see Hut 1982 for the dependence of the different factors on eccentricity). I don't consider this 
factor here, because the mass loss process, if it depends on the orbital phase, may lead to a change 
in eccentricity, which is not considered here either. The giant slows-down because of angular 
momentum loss via mass loss, at a rate given by 

4 = PRlcoM^, (5) 

where uj is the giant's angular velocity, is the mass loss rate in the wind, and /? is a parameter 
that depends on the mass loss geometry; /3 = 2/3 for spherical mass loss and /3 = 1 for equatorial 
mass loss. Both Jw and Mw are defined positively. In systems in which the evolution time is less 
than the circularization time, the angular velocity of the giant is given by equating the rate in 
equation (4) to that in equation (5), where I neglect the nagular momentum of the spining giant's 
envelope. This gives 



' Mg + AI, \ (a^ ^ (l-e^/^ l 
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where a;Kep is the Keplerian angular velocity on the equator of the giant, and /x = MgM2/{Mg + 
M2) is the reduced mass of the binary system. Substituting typical values and omitting terms of 

weak dependence on several factors in equation (3) (in any case, as indicated above, Tdr serves as 
approximation for the rate of angular momentum transfer in eccentric orbits) and taking M2 — 
Mg ~ M0, this ratio is 

— =.0.03/(e)(^— j j ( io-»Meyr- j ' 

where /(e) = /s(e)(l — e^)^/^ is a steeply increasing function of e, with /(O) ~ 1. I emphasize again 
that this relation holds only if the time spent by the mass losing star in the giant phase is shorter 
than the synchronization time Tgyn, only for mass loss via a wind (not for a Roche lobe overflow), and 
only if it gives uj < ujo, where Uo is the orbital angular velocity. The synchronization time is related 
to the circularization time by the expression Tgyn ^ (1 + M2/Mi){M2/Mi)~^{Ig/MgB?g){R/afTcir 
(accurate for e = 0), where Ig is the giant's moment of inertia. Using the typical values as in 
equation (3), I find equation (7) to hold for 3i?p ^ a ^ 20i?g, for a circular orbit; the separation 
can be larger for eccentric orbits. 

Crudely, equation (7) indicates that in all symbiotic systems with orbital period of Porb ^ 
100 yr, with larger periods for larger eccentricity, tidal interaction overcomes angular momentum 
loss, and the systems are synchronized, or close to being so. For systems with a ^ 5i?p, the situation 
is reversed during the superwind phase at the termination of the AGB, when the mass loss rate 
can be as large as ~ 10~^Mq yr"-*^, and the giant envelope loses synchronization and slows down. 
This is more relevant to the formation of some bipolar planetary nebulae than to SSs. 



2.3. Backfiowing material 



The basic process here is that the hot companion prevents, directly or indirectly, part of the 
mass blown by the giant companion to acquire the escape velocity from the binary system. That 
fraction of mass may acquire orbital angular momentum, and if it is accreted back by the giant, it 
spins-up the giant's envelope. This process should be worked out in detail via numerical simulations, 
e.g., similar to those of Mastrodemos k, Morris (1999), but with more processes included, e.g., locally 
enhanced mass loss rate. The conditions for a backflow to occur were examined in an earlier paper 
(Soker 2001). I found there that for a significant backflow to occur, there should be a slow dense 
flow, such that the relation between the total mass in the slow flow, Mj, and the solid angle it covers 
O, is given by Mj/(0/47r) ^ Q.IMq. The requirement for both high mass loss rate per unit solid 
angle and a very slow wind, such that it can be decelerated and flow back, probably requires close 
binary interaction. This process, therefore, requires that the companion be close and already some 
spin-up of the giant's envelope have occurred. Large magnetic cool spots (see next section) may 
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then lead to a very slow and a high mass loss rate above these spots. Because the escape velocity 
from the binary system is larger than that from the mass losing star, some of this material may 
be accreted back by the giant (or the companion). The gravitational interaction with the binary 
system may transfer some orbital angular momentum to the backflowing mass, e.g., in the accretion 
column formed behind the companion. The backflowing material, then, may possess large specific 
angular momentum and spin-up the giant. 

I suggest this process to account for the rotation of the SiO maser shell in R Aquarii (Mollis et 
al. 2001). Hollis et al. (2001) argue for a rotational period range of 8-34 yr, at a radius of ~ l — 2Rg, 
where Rg = 1.8 AU is the giant's radius. The orbital period of the R Aquarii binary system is 
^ 44 yr, with a semimajor axis of ~ 17 AU, and an eccentricity of ~ 0.8 (Hollis et al. 2001). The 
tidal synchronization time is very short for this system, hence, the angular velocity of the spinning 
giant should have been ~ 44 yr. Any mass lost by the giant should have a longer rotation period, 
and not as short as ~ 8—34 yr. The fast rotating maser spots are formed by backflowing material, as 
is indeed observed by Boboltz, Diamond, &; Kamball (1997). In particular, a high mass loss rate is 
expected during the periastron passage. A substantially enhanced mass loss rate during periastron 
passage may account for the high eccentricity of the system, despite the short circularization time 
(Sokcr 2000a). The backflowing mass, which is the source of the SiO maser, further spins-up the 
giant's envelope. Also, because of the fast rotation of the backflowing gas, there will be a shear as 
it reaches the stellar surface. Such a shear may further amplify the magnetic field near the equator 
(Soker 2000b). 



3. MAGNETIC ACTIVITY 

The three processes discussed in the previous section show that the giant companion in many 
SSs is likely to rotate much faster than isolated, or in wide binary systems, RGB and AGB stars. In 
the present section I discuss plausible implications of this (relatively) fast rotation to the magnetic 
activity of these cool giants. 



3.1. Global mass loss geometry 

A strong magnetic activity may influence the mass loss geometry, e.g., cause a higher or lower 
mass loss rate in the equatorial plane, and/or cause semi-periodic variation via a solar- like magnetic 
activity cycle. The magnetic field may influence the mass loss process via dynamical effects, i.e., 
the magnetic force due to pressure gradient and/or tension directly influences the outflowing mass 
(e.g.. Chevalier & Luo 1994; Pascoli 1997; Matt et al. 2000; Blackman et al. 2001; Garcia-Segura et 



- 7- 



al. 1999; Garci'a-Segura, Lopez, & Franco 2001; Gardiner & Prank 2001). As far as the shaping of 
circumstellar gas is concerned, my view is (Soker & Zoabi 2002; Soker 2002) that the global shaping 
mechanisms arc based on direct influence of a companion (e.g., Mastrodemos &; Morris 1999) or 
indirect effects of the magnetic field but not on direct dynamical effects of the giant's magnetic 
field. I find this to be the case in the SS system R Aquarii, from the following consideration. In the 
sun the magnetic activity determines the mass loss rate. Indeed, the average solar X-ray luminosity 
resulting from the magnetic activity is within an order of magnitude of the kinetic energy carried by 
the solar wind. In the SS R Aquarii, the X-ray luminosity of the central source is much lower than 
the rate of kinetic energy carried by the giant's wind. From the results of Kellogg, Pedelty, &: Lyon 
(2001), the central X-ray luminosity is L^^ 4 x 10^^ erg s~^. The rate of wind's kinetic energy 
is = 3 X 10^° erg s^^(Mw/lO~''M0 yr"^)(v„/10 km s^^)^, where Vyj is the wind speed. The 
mass loss rate from R Aquarii is somewhat smaller than the scaling above (see data summarized 
by Iben & Tutukov 1996), but the outflow speed is larger (Hollis et al. 1999). The low ratio of 
Lx/L^ < 0.01 in R Aquarii hints that the mass loss process from the cool giant is not determined 
directly by magnetic activity. Indirect effects of the magnetic fleld are very likely, though. In §3.3 
I speculate that some of the central X-emission of R Aquarii may result from magnetic activity on 
the giant's surface. 

An indirect effect can be the formation of magnetic cool spots, above which dust formation, 
hence mass loss rate, is enhanced (Soker 2000b and references therein). Globally, the magnetic 
energy is much below that carried by the wind, although in local spots the magnetic fleld can be 
strong. Magnetic cool spots can be concentrated in the equatorial plane as in the sun, leading to 
higher mass loss rate in the equatorial plane. It is also possible, probably in rapidly rotating stars, 
that the spots are concentrated in the polar regions, as in rapidly rotating MS stars (e.g., Schrijver 
&; Title 2001). To conclude, relatively strong magnetic activity may lead to axially symmetric mass 
loss process. 

If a magnetic activity cycle exists, then the mass loss rate may be semi-periodic (Soker 2000b). 
This means that a mass loss rate presently observed in a speciflc system may not represent the 
average mass loss rate. 



3.2. Stochastic mass loss 

One manifestation of magnetic cool spots is locally enhanced mass loss rate, i.e., in particular 
directions the mass loss rate is higher for a short time. A similar effect can be caused by spots 
formed by large convective elements (Schwarzschild 1975). This means, as mentioned above, that 
the surface of the cool giant can be asymmetric. Another implication is a stochastic mass loss rate in 
the direction of the hot companion, hence the accretion rate and geometry onto the hot companion 
may vary during short periods, i.e., weeks to years (much shorter than a possible magnetic cycle). 
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This stochastic accretion process may have the following observable manifestations. (1) Impulsive 
jets: A high mass flux in a specific direction toward the companion, but not exactly on the line to 
it, means that the flow has a high spcciflc angular momentum about the hot companion. If a disk 
was not present before, a temporary disk may be formed. Such a disk may blow two jets for a short 
time. (2) Density variation along the jets' axes: If a disk was present with two jets, the enhanced 
accretion may lead to a denser jet. The consequence is two continuous jets, but with variation in 
density along the jet axis. (3) A tilted and precessing disk: The enhanced mass flux may be away 
from the equatorial plane. In such a case, the angular momentum accreted may have a direction 
tilted with respect to the permanent accretion disk. This may tilt the disk and cause a precession 
for a short time. This may be observed as precessing jets, but for a limited time. 

Variation in the mass accretion rate may result from the pulsation of the giant as well, but to 
a lesser extent than I expect in a strongly magnetically active giant. 



3.3. X-ray emission 

Magnetically active MS stars have X-ray luminosity of up to ~ 10^° erg s~^ (Stelzer & 
Neuhauser 2001). The source of the X-ray luminosity is a hot corona and magnetic reconnec- 
tion events, e.g., flares. Some observations indicate that cool giants may also emit in the X-ray 
band (Hiinsch ct al. 1998; Schroder, Hiinsch & Schniitt 1998; Hiinsch 2001), resulting most likely 
from magnetic activity. The presence of magnetic flelds in AGB stars is inferred from maser po- 
larization (e.g., Kemball & Diamond 1997). It seems that even a very slow rotation is enough to 
sustain a turbulent dynamo activity in AGB stars, i.e., an o?u) dynamo (Soker h Zoabi 2002). 
The fast rotating cool giants in SSs may be among the most magnetically active cool giants and 
those in which strong X-ray emission is expected. There is no basic dynamo theory to predict the 
level of magnetic activity in AGB stars; crude estimates for AGB stars rotating with periods of 
~ 10^ — 10^ yr give Eb < 10^° — 10^^ erg s^^ (Soker & Harpaz 1992). The temperature of the X-ray 
emitting plasma is most likely ~ 10^ K, as in MS stars. Magnetic flares may lead to detectable 
X-ray variability on time scales of weeks to months. I therefore speculate that the X-ray central 
source in the SS R Aquarii, with Lj, ~ 4 x 10^^ erg s"^ and = 2 x 10^ K (Kellogg et al. 2001) 
result from a magnetic activity on the surface of the Mira variable. X-ray observations of other SSs 
is highly encourage. 



4. SUMMARY 

The main goal of the present paper is to point out that the cool giant companion in most 
symbiotic binary systems (SSs) are likely to be strongly (relative to other cool giants) magnetically 
active. 
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The observational aspects of the magnetic activity are: (1) Imposing a global axisymmetrical 
mass loss geometry from the cool giant, i.e., a higher mass loss rate in the equatorial or polar 
directions. (2) Stochastic mass loss process, i.e., variation of mass loss rate with time and location 
on the giant's surface. The hot companion itself may cause a much larger departure from spherical 
symmetry in the circumbinary matter. Therefore, the asymmetry due to magnetic activity, both 
global and local, should be looked for near the giant's surface, e.g., SiO maser spots. The stochastic 
mass loss process may also cause variations in the intensity and direction of jets blown by an 
accreting companion. (3) As in main sequence stars, the magnetic activity may lead to X-ray 
emission. The expected temperature is ~ 10^ K. The X-ray luminosity is hard to predict as there 
is no basic dynamo model for cool giants. A crude estimate is < lO^'^ — 10^^ erg s~^, in any 
case below the kinetic energy of the giant's wind. I suggest that this may contribute to the central 
X-ray source in R Aquarii, and hence I encourage more X-ray observations of SSs. 

In arguing for magnetic activity of the cool giants of SSs, I assumed, based on the behavior of 
main sequence stars, that magnetic activity and X-ray luminosity increase with rotation velocity. 
I then showed that the cool companions in SSs are likely to rotate much faster than isolated, or in 
wide binary systems, cool giants, because of three processes: (1) In SSs in which the hot companion 
is a WD, the presently giant accreted mass from the wind of the WD's progenitor during the AGB 
phase of this progenitor. The accreted mass has a large specific angular momentum, hence it spins- 
up the MS progenitor. (2) The giant is spun up via tidal interaction with the hot companion. (3) 
Some fraction of the mass blown by the giant may flow back to the giant after acquiring angular 
momentum from the binary system, e.g., gravitational interaction with the companion. 

This research was supported by the US-Israel Binational Science Foundation. 



-10- 



REFERENCES 

Belczynski, K., Mikolajewska, J., Munari, U., Ivison, R. J., & Priedjung, M. 2000, A&AS, 146, 407 

Blackman, E. G., Prank, A., Markiel, J. A., Thomas, J. H., & Van Horn, H. M. 2001, Nature, 409, 
485 

Boboltz, D. A., Diamond, P. J., & Kemball, A. J. 1997, 487, L147 
Chevalier, R. A., & Luo, D. 1994, ApJ, 421, 225 

Garcia-Segura, G., Langer, N., Rozyczka, M., &; Pranco, J. 1999, ApJ, 517, 767 
Garcia-Segura, G., Lopez, J. A., &; Pranco, J. 2001, ApJ, 560, 928 
Gardiner, T. A., & Frank, A 2001, ApJ, 557, 250 

Hollis, J. M., Bobolotz, D. A., Pedelty, J. A., White, S. A., & Forstcr, J. R. 2001, 559, L37 

Hollis, J. M., Vogel, S. N., Van Buren, D., Strong, J. P., Lyon, R. G., k, Dorband, J. E. 1999, ApJ, 
522, 297 

Hiinsch, M. 2001, in Astronomical Gesellschaft Abstarct Scr., Vol. 18, MS 07 10 

Hiinsch, M., Schmitt, J. H. M. M., Schroder, K.-R, & Zickgraf, P.-J. 1998, A&A, 330, 225 

Hut, P. 1982, A&A, 110, 37 

Iben, L Jr., & Tutukov, A. V. 1996, ApJS 105, 145 

Karovska, M., Hack, W., Raymond, J., & Guinan, E. 1997, ApJ, 482, L175 
Kellogg, E., Pedelty, J. A., & Lyon, R. G. 2001, ApJ, 563, L151 
Kemball, A. J., & Diamond, P. J., 1997, ApJ, 481, Llll 

Marengo, M., Karovska, M., Pazio, G. G., Hora, J. L., Hoffmann, W. P., Dayal, A. k, Deutsch, L. 
K. 2001, ApJ, 556, L47 

Mastrodemos, N., & Morris, M. 1999, ApJ, 523, 357(MM99). 

Matt, S., Balick, B., Winglee, R., & Goodson, A. 2000, ApJ, 545, 965 

Miranda, L. P., Gomez, Y., Anglada, G., & Torrelles, J. M. 2001, Nature, 414, 284 

Morris, M. 1987, PASP, 99, 1115 

Pascoli, G. 1997, ApJ, 489, 946 

Schrijver, C. J., & Title, A. M. 2001, ApJ, 551, 1099 
Schwarzschild, M. 1975, ApJ, 195, 137 

Schroder, K.-R, Hiinsch, M., & Schmitt, J. H. M. M. 1998, A&A, 335, 591 
Soker, N. 1998, MNRAS, 299, 1242 
Soker, N. 2000a, A&A, 357, 557 



- 11 - 



Soker, N. 2000b, ApJ, 540, 436 

Soker, N. 2001, MNRAS, 324, 699 

Soker, N. 2002, preprint ( ^stro-ph/0204157] ) 

Soker, N., & Harpaz, A. 1992, PASP, 104, 923 

Soker, N., & Zoabi, E. 2002, MNRAS, 329, 204 

Stelzer, B., k Neuhauser, R. 2001, A&A, 377, 538 

Stencel, R. E., & Sahade, J. 1980, ApJ, 238, 929 

Szymczak, M., Cohen, R. J., & Richards, A. M. S. 1999, MNRAS, 304, 877 

Verbunt, F., & Phinney, E. S. 1995, A&A, 296, 709 

Wang, Y.-M. 1981, A&A, 102, 36 

Zahn, J-P. 1977, A&A, 57, 383; erratum 67, 162 

Zahn, J-P. 1989, A&A, 220, 112 

Zijlstra, A. A., te Lintel Hekkert, P., Pottasch, S. R., Caswell, J. L., Ratag, M., & Habing, H. J. 
1989, A&A, 217, 157 



This manuscript was prepared with the A AS lATgX macros v3.0. 



